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CODED LABEL INFORMATION EXTRACTION METHOD 

The present invention relates to methods of 
interpreting coded information carriers and to methods 
5 of processing the response of such carriers to an 
interrogating magnetic field. 

In certain types of human and machine-readable 
information-bearing label, a set of elements is used to 
represent the information contained in the label. This 
10 representation may be made by varying the 

characteristics of the elements comprising the label, 
and also by the position in which the label elements 
are placed. Reading apparatus senses the 
characteristics and placement of' the elements in order 
15 to decode the information contained within the label. 

The elements in the label are sequentially scanned, to 
discover the presence of the constituent elements and 
to measure their characteristics and position and 
thence to decode "the information contained by the 
20 label. 

This application describes label decoding methods 
by which information represented on such a label can be 
decoded. Typical label embodiments are where the label 
is manufactured from a plurality of magnetically active 
25 elements supported on a substrate. Information is coded 
by controlling the relative positions. Independent 
control of the physical properties of the elements, 
such as shape may also be used. 

In alternative embodiments, information may be 
30 coded by controlling the relative angles of orientation 
of the easy axis of the magnetic elements and in some 
cases element shape some other physical property can be 
independently controlled to further encode the 
information carrier. 
35 in certain types of label decoding system a label 

consisting of a set of elements is scanned, that is the 
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elements are activated in turn by applying a signal to 
which the scanned element responds, and to which the 
reading apparatus is sensitive. Either or both of the 
means of generating the scanning signal, and the means 
5 of generating the received signal in the reading 

apparatus select a single elements in turn so that the 
signal from a particular element may be resolved from 
that of the other elements. 

According to one aspect of the present invention 

10 there is provided a method of decoding an information 
carrier comprising a plurality of magnetic elements 
supported by, or incorporated in, a substrate, wherein 
the relative positions and/or physical dimensions of 
said magnetic elements represents the information 

15 encoded by said information carrier, the method 

comprising the steps of: 

i) applying an interrogation signal to each of the 
elements of said information carrier; 

ii) detecting the response of said magnetic 
20 elements to said interrogation signal? 

iii) processing the response detected in ii) so as 
to determine an estimate of a mathematical function 
which defines the relative velocity between the 
interrogating signal and each of the magnetic elements, 

25 from which the relative positions and/or the relative 

physical dimensions can be determined. 

Conveniently, detection of the magnetic response 
of said magnetic elements may comprise observation of 
harmonics which are generated by the magnetic element 

30 from an applied AC field as its magnetic state is 
altered by the interrogating magnetic field . The 
presence of a magnetic element, and details of the 
physical dimensions of the magnetic element, can 
advantageously be identified in a plot of the detected 

35 amplitude of the harmonic signal. The presence of a set 

of signals in the otherwise quiescent detected signal 
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indicates the issuance of harmonic signal by a magnetic 
element as its magnetic properties are active. The 
means of scanning may be automatic by electronically 
scanning the interrogating signal and the receive 
5 sensitivity along the label or by an automatic 

mechanism causing the label elements to be interrogated 
in turn. Alternatively the scanning may be manually 
operated where the label is moved past the reading 
apparatus by hand, or alternatively the reading 
10 apparatus or part thereof may be moved past the label 
by hand. 

Information may be advantageously encoded into the 
label by varying either the length of the material 
elements, that is the size of the element normal to the 

15 longitudinal axis of the label, or the width of the 

elements. By maintaining the element width constant, 
variations in element length predominantly cause a 
change in amplitude of the received signal which can be 
detected in order to recognise additional coded 

20 information. 

For a better understanding of the present 
invention and to show how the same maybe carried into 
effect, reference will now be made by way of example to 
25 the accompanying drawings in which: 

Figure la shows a coded information carrier 
according to the present invention; 

Figure lb shows the variation in amplitude of the 
second harmonic frequency as the information carrier is 
30 subjected to an interrogating magnetic field; 

Figure 1c shows the receive signal from a single 
magnetic element; 

Figure Id shows the received fundamental signal 
from three magnetic elements; 
35 Figure 2a shows the relative variation of the 

received second harmonic signal with element length; 
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Figure 2b shows the relative variation of received 
fundamental signal with element length; 

Figure 3a show a coded information carrier 
according to the present invention, wherein the 
5 magnetic element are length coded; 

Figure 3b shows the receive signal from an 
information carrier having length coding; 

Figure 4a shows the positions of the half peaks of 
the detected signal shown in Figure 3b; 
10 Figure 4b shows the magnitude of the double peaks 

in amplitude units ; 

Figure 4c shows a plot of the double peak 
separations versus double peak position; 

Figure 4d shows a plot of the velocity estimates 
15 of selected points versus time; 

Figure 4e shows the resultant fit of gap number to 
integer grid. 

Fig la shows a plan view of a label 1 comprising a 
20 plurality of magnetic elements 2. The elements are made 
out of high permeability, low coercivity magnetic thin 
film material such as Atalante manufactured by 
1ST (Belgium) . In this particular example the elements 
are substantially equal in size and have dimensions of 
25 lmmxl.Smm. The easy axis of the material is arranged so 

as to be parallel with the longitudinal axis of the 
label* 

In order to interrogate the label 1, a reading 
antenna is passed by the label subjecting each element 

30 in turn to a region of null magnetic field formed by a 
permanent magnet arrangement. That is a magnetic field 
comprising regions of high saturating magnetic field 
contiguous with regions of zero magnetic field. A 
superimposed high frequency field, in this case 6.348 

35 kHz, causes the material element to switch when it is 

present in the magnetic null. In other words the 
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element is driven over its B-H loop as it passes 
through the zero magnetic field region, A receive coil 
is used to detect the harmonic signals resulting from 
the element switching, and the received harmonic signal 
5 is mixed to baseband and filtered prior to 

digitisation. 

Figure lb shows a typical signal received from 
such a label, wherein the read head is passed in 
proximity with the information carrier and is arranged 

10 to detect the second harmonic of the element switching 
signals. The signal from a single element is shown in 
Fig 1c which shows two half peaks 4 and 5 of opposite 
sign with a zero crossing. It can be seen that the 
label signal approximates a superposition of single 

15 element responses at each position in which there is an 
element . 

The signal is sampled with a fixed sample rate, 
and converted to a digital representation in an 
analogue to digital converter. The sample rate must be 
20 fast enough to capture the signal information in 

accordance with the Nyquist criterion, preferably with 
anti-aliasing filtering prior to the sampling 
operation . 

The sample number of the centre of each element 
25 response depends upon the time at which the element 
passed the read head. 

In addition, the width of the element response 
depends upon the speed at which the element passed the 
read head. Information about the dynamics of the 
30 relative motion of the tag and the read head can 

therefore be inferred from measuring the width of the 
signal. This can be achieved by measuring the width of 
the positive going and negative going peak associated 
with the element. It can also be achieved by measuring 
35 the separation between the positive and negative going 
peaks. This speed estimate is fundamental to estimates 
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of the dynamics of motion and allows the information in 
the label to be decoded reliably without compromising 
the available code space by adding regular features in 
the label. 

5 In alternative implementations of the reading 

apparatus, the scanning motion of the magnetic null 
relative to the label may be achieved by electronically 
scanning using a set of excitation coils that are 
driven in a phased manner so as to cause the magnetic 

10 null to propagate. in space- The null may be scanned at 
a speed which allows the switching of the label 
elements to be detected directly from the receive coils 
without needing to superimpose a separate high 
frequency excitation field. In this case the material 

15 element response coupled into the receiver is a single 

peak. 

Fig Id shows 3 adjacent peaks 6, 7 and 8 of the 
fundamental signal from 3 magnetic elements recorded 
from a reading system employing electronic scanning. In 

20 this case the width of the response from each magnetic 
element can be used to provide an estimate of the 
velocity of the scanning null as it passed over the 
element being scanned. This information can be used to 
estimate the detailed dynamics of motion and the errors 

25 arising from the system and also from misalignment of 
the label with respect to the nominal read axis of the 
read head. 

Information can be encoded into the label by 
varying the length of the material elements, that is 
30 the size of the element normal to the longitudinal axis 
of the label, or the width, that is the size of the 
element along the longitudinal axis. 

By maintaining the element width constant, 
variations in element length predominantly cause a 
35 change in amplitude of the received signal as 

illustrated in Fig 2a for the mechanically scanned 
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second harmonic system, and in Fig 2b for an 
electronically scanned system detecting the fundamental 
signal received from the label. 

A set of linear test labels with 3 elements was 
5 manufactured with the element positions at Omm, 1.5mm, 
3mm along the longitudinal axis of the label. The 
element widths were 1mm and the length of the centre 
element varied between linm and 2mm. For each label, in 
Fig 2a, the peak position, the half peak width, the 

10 half peak separation and the amplitude difference 

between the half peaks is advantageously measured. The 
relative variation of these parameters of the central 
test element compared with the outer elements, which 
were maintained at constant length, may be plotted as 

15 shown in figure 2b. The measured central element 

amplitude (vampr) , relative to the outer elements, 
shows a clear, near linear relationship to central 
element length. The measured central element relative 
position (vposr) , width (vwidr) and separation (vsepr) 

20 parameters do not vary significantly with central 
element length, over and above the effects of 
experimental errors and imperfections in the 
manufacture of the test labels - 

Similarly Fig 2b shows the relative variation of 

25 measured central element amplitude, position and width 
with element length when the test labels were measured 
in an electronically scanning reading system and shows 
a similar near linear relationship to the mechanically 
scanned reading method receiving second harmonic 

30 signals. 

A label in which element length has been used to 
encode additional information is shown in Fig 3a , 
which shows a label 10 in which two element [with 
elements of one of two] lengths are used to modulate 

35 the elements in the code. The corresponding received 
signal is shown in Fig 3b in which a difference in 
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received amplitude from elements of different lengths 
can be observed. A typical design goal is to maximise 
the range at which a label may be read. A compromise is 
needed to maximise reading range for a 2 -length label, 
at the same time as maximising the ability to 
discriminate the lengths of individual elements. In 
order to discriminate the length of label elements 
reliably, certain of these elements may be constrained 
to be at a known length to act as a reference. In the 
example illustrated in 3a, the 2 end elements, 12 and 
13 , of the label are maintained at maximum length to 
act as a reference. A design rule may also constrain a 
certain number of the inner elements to also be at 
maximum length. The reference elements may be at any a 
priori fixed position, or at a position determined by 
an a priori rule. 

The half peak width, or half peak separation can 
be measured and, in the case of a reader detecting the 
second harmonic of the label signal, can be used to 
form an estimate of the velocity of the interrogation 
field, at the time when a particular element was 
scanned, Alternatively, in the case of a reader 
detecting the fundamental of the label signal, an 
estimate of the peak width can be used to form an 
estimate of the velocity of the scanning, at the time 
when the element was scanned. These estimates are 
typically K/w_est where K is some constant, and w_est 
is the estimate of width or half peak separation, as 
appropriate. The value of K typically depends on 
parameters of the reader, the element shape and the 
position of the element relative to the read antenna. 
Since K varies with element properties and to 
facilitate the estimation of the scanning dynamics, K 
may only be estimated to for elements of the same type 
based on the amplitude of the element signals. 

The first step in processing a received signal 
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such as shown in Fig lb or Fig 3b is to estimate the 
parameters of the peaks contained within the signal. 
This is achieved with the following steps: 

5 1 The signal is delimited, that is the continuous 

received signal is broken up into a contiguous set of 
samples containing the received signal from one 
label, or a set of peaks that is likely to contain 
enough information to decode the label, in the case 
10 of non- registered repeating codes, 

2 The signals from this delimited period cure sorted 
into order of ascending magnitude. An indexed sort 
function, for example using the Quicksort method, 
such as index 0 (Numerical Recipes in C, second 

15 edition, ISBN 0 521 43108 5) is appropriate so that 

the original location of the point of any rank can be 
quickly found- The median value is found from the 
sorted data and provides a good estimate of the 
quiescent signal level. The median is subtracted from 

20 all signal points so that the quiescent level is 

close to zero. 

3 The positive half peaks are found by working down 
from the largest signal in magnitude in the sorted 
data, testing the signal points as follows. If the 

25 test signal point is a peak, that is its adjacent 

signals are lower in magnitude it is estimated as a 
positive half peak, and its amplitude is interpolated 
by interpolation of the peak and its neighbouring 
points . 

30 4 From each positive half peak found, the signal is 

scanned forward and backwards to find the points at 
which it falls below a defined magnitude relative to 
the peak, this is interpolated to a fraction of a 
sample period A typical magnitude value is 0.8 of the 

35 peak value found. The half peak width is estimated as 
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the difference in position between these points . The 
half peak position is estimated as the average of 
these two points. The half peak is rejected if it 
overlaps with any other pre-existing peak v in the 
5 list, that is its width about its position interferes 

the width about any other peak position in the list. 
If it is not rejected, the half peak is appended to a 
list of half peaks. 

5 The peak finding process stops when the magnitude of 
10 the next highest point falls below an absolute 

threshold, or when it falls below a defined fraction 
of the largest magnitude signal, for example 0.3 of 
the largest magnitude signal, or if a limit for the 
number of peaks is reached. 

15 6 The negative half peaks are found by working up in 

magnitude from the lowest signal in magnitude in the 
sorted data, testing signal points of increasing 
magnitude. If the test signal point is a negative 
peak, that is its adjacent signals are higher in 

20 magnitude it is estimated as a negative half peak, 

and its amplitude is interpolated by interpolation of 
the peak and its neighbouring points. 

7 For 'each negative half peak, the signal is scanned 
forward and backwards to find the points at which it 

25 increases above a defined magnitude relative to the 

negative peak, this is interpolated to a fraction of 
a sample period. The half peak width is estimated as 
the difference in position between these points. The 
half peak position is estimated as the average of 

30 these two points. The half peak is rejected if it 

overlaps with any other pre-existing peak in the 
list, that is its width about its position interferes 
the width about any other peak position in the list. 
If it is not rejected by this test it is then 

35 appended to a list of half peaks. 
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8 The negative peak finding process stops when the 
magnitude of the next highest point falls above an 
absolute threshold, or when it falls above a defined 
fraction of the largest magnitude signal, or if a 

5 limit for the number of peaks is reached. 

9 The resultant list of positive and negative half 
peaks is then sorted into range order. This process 
will interleave the consecutive half peaks, 

10 Merging range consecutive pairs of half peaks in 
10 the range sorted half peak list forms a list of 

double peaks. The amplitude of the merged double peak 
is the difference in amplitude of two half peaks 
being merged. The position of the merged double peak 
is the average of the two half peak positions. The 
15 width of the merged double peak is the average of two 

half -peak widths. The separation of the merged double 
peak is the difference between the two half peak 
positions. 

20 The results of this half peak finding process are 

illustrated in Fig 4a which shows the half peaks found 
when the signal shown in Fig 3b is processed, and shows 
the output of step 9 above. These peaks are shown with 
a cross at the located position and estimated 

25 amplitude. Fig 4b shows the magnitude of the double 

peaks in amplitude units. In this case the amplitudes 
are negative since the first half peak is a negative 
going peak. The difference in peak amplitude between 
the full length and reduced length elements in the 

30 label can be clearly seen. Fig 4c shows the double peak 

separations versus double peak position showing the 
variation across the scanning of the label. 

In the case of a fundamental signal where there is 
just one peak per element, such as that shown in Fig 

35 id, a simpler process is may be used, missing steps 6 
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to 10 above since there is no need to find and merge 
double peaks. 

There are N label elements positions whose centres 
are located at X t . . X» along the longitudinal axis of 
5 the label/ and whose position we wish to estimate. It 

is assumed that Xj is at position 0.0 and that the final 
element, X„ is at the known position of corresponding to 
the length of the label L. 

The peak list as described above contains 

10 estimates of the sample number, and hence times T x at 
which the reader detected a double peak. For each 
element i, the time T*. is simply the position in 
samples where the peak was measured multiplied by the 
sample period, which is the reciprocal of the sampling 

15 rate. 

By using the estimated double peak width or double 
peak separation for each of the N peaks denoted W x a 
velocity estimate for the i th peak is VE^lM at time 

20 The order M velocity function of time is modelled 

by a polynomial f unction. M is chosen to allow 
sufficient degrees of freedom to model the dynamics of 
label motion sufficiently accurately. It can range from 
1 (constant acceleration model) upwards . 

25 

VF(t,M) =a 0 +a a t+a 2 t 2 +a 3 t 3 +....a M t M 

where {a} are a set of polynomial coefficients to be 
determined. 

30 

The squared error for velocity estimate VE t at time T d 
is E x 2 = (VEi-VF (T^M) ) 2 

The total sum of the squared errors is ES(M)=E! 2 + 

35 E 2 2 + + E H 2 which is minimised by varying the M+l 

parameters a 0 ,a a a„ This is achieved by the method of 
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least squares to find the optimum parameters in the 
least squares sense. 

Define the matrix X to contain the T* values raised 
to a given power. For example column 0 row i of X 
5 contains and column M row 4 of X contains . 

Define the vector V to contain the velocity estimates 
VEi . For example the element in row 4 of V contains VE 4 

Some linear combination of the columns in X will 
be closest to the vector V in the least -squares sense- 

10 The coefficients of this linear combination will be the 
coefficients of the best-fitting polynomial. If A is 
the vector of length M+l that holds the coefficients, 
the problem is to find the value for A that gives the 
best "near -solution" to the system of equations XA= V. 

15 If the columns of X are independent, the best choice 

for A is given by: 

A=<x T xrMx T v) 

20 Where " T " represents matrix transpose and n " l w 

represents matrix inversion. The coefficient a x is then 
the element from row i of A. This may be efficiently 
calculated by using methods such as LU decomposition, 
for examples the functions ludcmpO and lubksbO 

25 (Numerical Recipes in C, second edition, ISBN 0 521 

43108 5) 

This process is illustrated in Fig 4d where the 
velocity estimates W a that were computed from the double 
peak separation data of Fig 4c are shown, those points 

30 marked "full data" contain the velocity estimates from 
all the peaks. In this example case, only the estimates 
from peaks with high amplitude, see Fig 4b, are 
selected for velocity fitting. These selected points 
are marked "Data points used in fit" . If the label 

35 known a priori to have elements of the same size then 

all elements would have been used. A least squares fit 
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of order 2 was used and the continuous line shows the 
estimated velocity function as a function of sample 
number. A least squares fit of order 1, 2 or 3 is 
typically sufficient to model the velocity function 
5 with sufficient accuracy. 

Error metrics are computed to allow the decoding 
algorithms to assess the quality of the read. A useful 
normalised velocity error metric, denoted evel, to 
quantify how well the velocity function fitted the data 
10 observed is defined. It is the sum of squares of the 
normalised error defined as : 

evel=sum( E^/VE*) for i=l..N 

15 The vector A could also be estimated by an 

iterative minimisation procedure such as Powell's 
method* 

The above procedure can be repeated for varying 
model order M on the same data in order to find the 
20 best approximation. Preferably this is the solution 
with the smallest number of degrees of freedom, that 
fits the data read by the reader with the lowest error 
metrics. 

25 Analytical interrogation and normalisation of tag 
length 

The velocity fit function is 

VF{t,M)=a 0 +a 1 t+a 2 t 2 +a 3 t 3 +.-..+a M t M This may be analytically 
30 integrated to find the positional fit function 

P(t,M) =Q+a 0 t+ (l/2)a!t 2 + (1/3) a 2 t 3 +.... + (1/ (M+l) )a„t Mn 
where Q is an arbitrary constant of integration, set to 
zero. 

35 

We know that the position of label element 1 is at 0.0. 
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We also know that the position of element N, the last 
element scanned is L, the length of the label. The 
final estimates of element position are: 

XEi-L . (P<X ± , M)-P(X a ,M)> / (PCX* ^J-PtX^M)) 

Binning of data 

There are N positional estimates XEi to XE„ from 
which N-l gaps may be estimated • 

Typically the gaps have the minimum gap 
subtracted, and are divided by the minimum gap, to give 
a real number representing the number of gap 
increments. The nearest integer Is then taken to 
represent the gap for subsequent decoding. This can be 
expressed as: 

The N-l normalised gaps are calculated by 
R 3 = (XE^-XE^MinimumJSap) /Gap_Increment for j =1„.N- 

1 

The normalised gaps are then converted to their nearest 
integer representaton 

G 3 = int(R>+0.5) for j=l.-N-l, where int(x) means the 

truncated integer portion of x. 

In this way, the gaps between each adjacent pair of 
magnetic elements can be defined by A + mG, wherein A 
is a first fixed value; m is an integer (which may be 
zero) ; and G is a second fixed value 

One useful error metric, denoted emax, is defined as 
the maximum magnitude of the difference between G 5 and 
R 3 .Another useful error metric, denoted eav, is the 
average magnitude of the difference between G 3 and R,. 
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Fig 4e shows the gap sequence data Ri computed from the 
data in Fig 4a-d, and shows the close correspondence or 
"snapping" to an integer grid, giving a high confidence 
that the gaps have been correctly measured. The gap 
5 sequence in this case is 3,2,2,1,7,0,0,2,1 . 

Enumeration algorithms 

Depending on the coding scheme used, the gap sequence 
is enumerated to a unique decoded value. Various codes 
exist that have different properties for example those 

10 codes allowing read direction to be determined, or 
those that allow non- registered operation, that is 
allowing the gap sequence to start at any position, for 
example a portion of a repeating gap sequence. The code 
may be determined a priori, or diagnosed by examining 

15 the read signals for the number of material elements 
present, or by estimating the sum of the gap numbers. 
It is preferable to determine this beforehand in order 
to minimise the chance of misinterpreting a given set 
of read data. 

20 

For certain codes, the number of elements in the label 
may be one of a range of possible values, with fixed 
label length in order to increase the number of 
possible codes. 

25 

The direction finding codes are preferably used where 
amplitude coding is used to allow the amplitude 
information to be associated with the correct element. 

Amplitude threshold estimation of elements 

30 Where it is desired to detect the size of the received 
element signal, for example when the elements may be of 
two or more sizes, reference elements may be placed in 
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the set of elements used to make the label. In one 
embodiment, the end elements are constrained to be the 
largest element size to act as a reference. Other 
embodiments may place reference elements that have some 
5 unique characteristic, or may be identified by an a 

priori rule. 

The amplitudes are measured and used to compute 
normalised amplitudes. This is because the reference 
10 elements may be significantly different in magnitude. 
This may be due to the scanning process varying the 
distance or angle of the read head from the label 
during the read, or other factors. 

15 The estimates of element position XE A may be used for 

linear interpolation of the reference amplitude along 
the label. The reference amplitude predicted at the 
estimated position of the element divides the measured 
element amplitudes to normalise the amplitudes. The 

20 normalised amplitudes should be nominally 1.0 for an 
element the same size as the reference element. 

Threshold levels are set to provide discrimination 
between the possible element amplitudes and to classify 

25 the amplitude of the element into one of an alphabet of 

values. For example if there are two values then the 
amplitude is either high or low. For example a half 
length element may have an amplitude, on average, half 
the reference amplitude ie 0.5. In this case, the 

30 threshold level should be set to approximately 0.75 in 
order to determine if the element was full length or 
half length. 

Preferably a code is used which enables the direction 
35 of reading of the label to be determined by the 

sequence of gaps that resulted. In this case the 
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elements can be uniquely identified and it is known 
which element of the label to associate with the 
decoded amplitude information. 

5 In the case of the label data shown in Fig 4b, denoting 
the amplitudes of the signals from the elements as D H" 
for high amplitude and "L" for low amplitude, the 
elements are categorised as H,L,H,L,H,Ii,H,L,L,H. A 
larger alphabet of elements could be used by using a 
10 label with more possible shapes, for example- The 

categorisations are then converted to decoded data. 

A typical scheme allows the inner elements only to 
contain data, ignoring the end elements, which are 

15 references, and maps H to a logical 0 and L to a 

logical 1 so that the binary data word from the example 
data would be 10101011. In this scheme up to 256 
amplitude codes can be represented from 8 inner 
elements. Preferably the codes used are constrained to 

20 ensure that a certain number of elements in the label, 
for example 4 are at full length to facilitate velocity 
estimation. If this constraint is used, the number of 
amplitude codes is reduced by 37 to 219 in this case. 

25 An amplitude error metric eamp is defined which is the 
normalised amplitude difference between the largest 
half length element diagnosed to the smallest full 
length element diagnosed. If there is no half-length 
element diagnosed, eamp is defined as the smallest full 

30 length element amplitude. In a label containing 

different sized elements, a large value for this metric 
suggests a clear discrimination in amplitude between 
elements of different types. Smaller values indicate 
increasing risk that an element may be categorised as 

35 the wrong amplitude. 
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Error Metrics 

The error metrics defined: evel, emax, eav and eamp can 
be used to determine the read quality, evel, emax and 
eav are preferably small and eamp high. Independent 
thresholds for these quantities are independently 
established. These thresholds are based on measured 
values for many labels in many read operations over the 
range of reading ranges and scanning dynamics expected 
in an application. This is done to minimise the 
probability of a wrong code at the same time as 
maintaining a high probability of a correct read. 

During the read decoding process, if the computed error 
metrics are higher than the threshold value in the case 
of evel, emax and eav then the read is rejected* If the 
error metric is lower than the threshold in the case of 
eamp then the read is rejected. 
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CLAIMS 

1. A method of decoding an information carrier 
comprising a plurality of magnetic elements supported 
5 by, or incorporated in, a substrate, wherein the 

elements are of substantially the same width and 
wherein the relative positions of said magnetic 
elements represents the information encoded by the 
information carrier , the method comprising the steps 
10 of: 

i) applying an interrogation signal to the 
information carrier so that each of the magnetic 
elements are subjected, in turn, to the interrogation 
signal; 

15 ii) detecting the response signal of each of the 

magnetic elements to said interrogation signal ; 

iii) processing each response signal detected in 

ii) so as to determine an estimate of the relative 
velocity between the interrogation signal and each of 

20 the magnetic elements; and 

iv) estimating the dynamics of motion between the 
interrogation signal and the information carrier, by 
consideration of the processed signals obtained in 

iii) , so as to determine the relative positions of the 
25 magnetic elements. 

2 A method according to claim 1, wherein said 
magnetic elements comprise high permeability, low 
coercivity material having an easy axis of 
30 magnetisation. 

3. A method according to claim 1 or 2, wherein said 
interrogation signal comprises regions of high 
saturating magnetic field contiguous with regions of 
35 zero magnetic field. 

J 

i 
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4. A method according to any preceding claim, wherein 
the frequency of the detected response generated by 
said magnetic elements is a harmonic of the 
interrogation signal. 

5. A method according to any preceding claim, wherein 
the frequency of the detected response generated by 
said magnetic elements is the second harmonic of the 
interrogation signal. 

6. A method according to any claim 4 or 5, wherein a 
graphical representation of the response signal is 
obtained by plotting the amplitude of the detected 
signal against time, wherein the presence of a magnetic 
element is indicated by two half peaks of opposite 
polarity in the detected response signal. 

7. A method according to claim 6, wherein the half 
peaks in the detected response signal are identified by 
the following steps: 

i) obtaining a contiguous set of sample points 
from the detected response signal; 

ii) sorting the samples points into an order of 
ascending magnitude; 

iii) finding the median value from the sorted 
samples and subtracting this value from each of the 
sample points; 

iv) testing the sample points in order to identify 
the peaks, wherein a positive polarity half peak is 
identified by the presence of lower magnitude signals 
either side of the signal point being tested, and a 
negative polarity half peak is identified by the 
presence of higher magnitude signal either side of the 
signal point being tested. 

8. A method according to claim 7, wherein the relative 
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velocity between the interrogation signal and each of 
the magnetic elements is determined by: 

i ) determining the two positions either side of 
the peak at which the amplitude of the detected 
5 response signal is a predetermined fraction of the peak 
amplitude value; and 

vi) determining the distance between these two 
positions to give the half peak width. 

10 9. A method according to claim 6 or 7, wherein an 

estimate of the relative velocity between the 
interrogation signal and the information carrier during 
its interaction with each of the magnetic elements is 
determined by measuring the separation between the two 

15 half peaks of opposite polarity. 

10. A method according to any one of claims 1,2 or 3, 
wherein the frequency of the detected response 
generated by said magnetic elements is the fundamental 

20 frequency of the interrogation signal. 

11. A method according to claim 10, wherein a 
graphical representation of the response signal is 
obtained by plotting the amplitude of the detected 

25 signal against time, wherein the presence of a magnetic 

element is indicated by a peak in the detected response 
signal . 

12. A method according to claim 11, wherein the peaks 
30 in the detected response signal are identified by the 

following steps: 

i) obtaining a contiguous set of sample points 
from the detected response signal; 

ii) sorting the samples points into an order of 
35 ascending magnitude; 

iii) . finding the median value from the sorted 



AMENDED SHEET 



-23- 

samples and subtracting this value from each of the 
sample points; 

iv) testing the sample points in order to identify 
the peaks, wherein a peak is identified by the presence 
5 of lower magnitude signals either side of the signal 

point being tested* 

13. A method according to claim 12, wherein an 
estimate of the relative velocity between the 

10 interrogation signal and each of the magnetic elements 

is determined by: 

i ) determining the two positions either side of 
the peak at which the amplitude of the detected 
response signal is a predetermined fraction of the peak 
15 amplitude value; and 

ii) determining the distance between these two 
positions to give the peak width. 

14. A method according to any one of claims 6 to 13, 
20 wherein the velocity of the interrogation signal, with 

respect to each magnetic element, is determined from 
the values of each of the peak or half peak widths. 

15. A method according to any preceding claim, 
25 wherein the velocity function of the interrogation 

signal is modelled by a polynomial function, the 
coefficients of which are determined by the estimates 
of the velocity of the interrogation signal at each 
magnetic element. 

30 

16. A method according to claim 15, wherein said 
mathematical function is integrated with respect to 
time so as to determine the size of the gaps between 
each of the magnetic elements. 

35 

17. A method according to any preceding claim, 
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wherein the size of the measured gaps between each 
adjacent pair of magnetic elements are modelled by A + 
mG, wherein A is a first fixed value; m is an integer 
(which may be zero} ; and G is a second fixed value . 

18. A method according to any preceding claim, 
further comprising the step of decoding a variation in 
length of the magnetic elements, wherein the amplitude 
of the detected signal response is directly related to 
the length of the element. 

19. A method according to claim 18, wherein the length 
of at least one of the magnetic elements is or are 
known, thereby acting as a reference element, with 
respect to which the lengths of the other magnetic 
elements can be determined. 

20. A method according to any one of claims 15 to 18, 
further comprising the step of determining an error 
metric of the velocity function. 

21. A method according to any one of claim 17, further 
comprising the step of determining an error metric of 
the measured gaps between the magnetic elements 
compared with the gap widths as modelled by A+mG where 
m is a decoded integer from each measured gap. 

22. A method according to any one of claims 15 to 18, 
further comprising the step of determining an error 
metric of the amplitude of the magnetic response. 

23. A method according to any one of claims 20, 21 or 
22, wherein the decoded information is rejected if the 
error metric is found to be higher than a threshold 
level. 
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